Ultrafast resonant optical scattering from single gold nanorods: 
Large nonlinearities and plasmon saturation 
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We measure nonlinear optical scattering from individual Au nanorods excited by ultrafast laser 
pulses on resonance with their longitudinal plasmon mode. Isolating single rods removes inhomo- 
geneous broadening and allows the measurement of a large nonlinearity, much greater than that of 
nanorod ensembles. Surprisingly, the ultrafast nonlinearity can be attributed entirely to heating of 
conduction electrons and does not exhibit any response associated with coherent plasmon oscillation. 
This indicates a previously unobserved damping of strongly driven plasmons. 

PACS numbers: 78.67.Bf,78.47.+p 
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I. INTRODUCTION 

Conventional photonic devices are restricted by the 
diffraction limit to be larger than half the optical wave- 
length, limiting the possibilities for miniaturization. 1 
One way to overcome this limit is to couple light to ma- 
terial excitations. Of particular interest are collective 
electron oscillations in metal nanostructures, termed sur- 
face plasmonsi^ Plasmons allow for nanoscale dereal- 
ization, transport of electromagnetic energy, and large 
local field enhancements. Preliminary steps have been 
made, for example, towards using surface plasmons in 
nanoparticles to construct sub- wavelength waveguides4£ 
Actively controlling light propagation in such structures 
will require an understanding of the ultrafast nonlinear 
response of the individual elements. Such an understand- 
ing is also crucial for the treatment of effects such as 
surface-enhanced Raman scattering(£*£ since it may limit 
the magnitude of local electromagnetic fields that can be 
achieved in real structures. 

Previous measurements of metal-nanoparticle non- 
linearities have generally involved excitation at fre- 
quencies away from the plasmon resonance, and thus 
probe incoherent effects related to the heating of 
electrons iSi2iifliiLiSiiii4 Exciting and probing nanopar- 
ticles on resonance with their plasmon frequencies can 
reveal nonlinearities associated with the coherent oscilla- 
tion of the plasmons themselves. Unfortunately, the op- 
tical response of the ensemble is broadened by the inho- 
mogeneous distribution of particle sizes and shapes. The 
majority of particles are off resonance with the exciting 
laser, and thus have nonlinear responses much smaller 
than or even opposite in sign to the resonant particles. 
This leads to an overall measured effect that is greatly 
reduced and whose dynamics are obscured. By contrast, 
isolating single particles allows the quantitative measure- 
ment of inherent propertiesi^iS^i 

We report the first measurements of resonant nonlin- 
earities of surface plasmons in single metal nanoparti- 
cles, specifically Au nanorodsiiSiiS We measure a nonlin- 
ear scattering cross-section that is much larger than that 



obtained from ensembles of nanorods. Surprisingly, the 
measured effect can be explained entirely as the result 
of heating of conduction electrons, with no measurably 
nonlinearity directly associated with coherent plasmon 
oscillation. This indicates that the strong optical driv- 
ing fields induce a novel damping and saturation of the 
plasmonic response. 



II. EXPERIMENTAL 

The Au nanorods we study are chemically synthe- 
sized using a seed-mediated growth methodi 20 ' 21 Un- 
der the proper growth conditions, this process produces 
single-crystal rods with smooth surfaces, controllable as- 
pect ratios, and > 95% yields A transmission-electron- 
microscope image of a typical nanorod is shown in Fig. 
^c). The rods exhibit a strong longitudinal plasmon res- 
onance, whose frequency is determined by the aspect ra- 
tio of the rods. 23 Damping due to interband transitions is 
reduced by selecting particles with a plasmon resonance 
near 1.55 eV^i which matches the TkSapphirc laser light 
used to excite and probe the rods. 

The gold- nanorod solution is prepared as follows. 
Spherical seed particles are produced by mixing 0.25 
mL of 10 mM HAuCl 4 solution with 10 mM of 0.1 
M CTAB (cetyltrimethylammoniumbromide) solution at 
room temperature, and then quickly injecting 0.6 mL of 
freshly prepared 10 mM NaHBH4 solution under vigor- 
ous stirring. In order to grow the seeds into rods, 50 mL 
of 0.1 M CTAB solution is prepared and maintained at 
28 °C. To this solution, 2.5 mL of 10 mM HAuCl 4 , and 
0.5 mL of 10 mM AgN0 3 are added. 1.0 mL of 1.0 M 
HC1 is also added, to maintain the stability of the final 
produce. Au(III) is reduced to Au(I) by injecting 0.4 
mL of 0.1 M ascorbic acid. Finally, 0.12 mL of the gold- 
seed solution is added to begin the nanorod growth, and 
growth is completed overnight under steady stirring. No 
further size selection is performed. 

The sample consists of sparsely dispersed and iso- 
lated rods, bound to a glass coverslip. The coverslip is 
first cleaned for 10 minutes with an equal mixture of 
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FIG. 1: (Color online) (a) Extinction of an ensemble of Au nanorods in aqueous solution (dashed blue line); scattering 
spectrum from a single nanorod on a glass surface (circles); calculated scattering spectrum for a single rod (solid black line); 
and measured spectrum of the laser used to excite the rods (solid red line), (b) Intensity of laser light scattered off a single rod 
as the incident laser polarization is varied (squares), and sinusoidal fit (red line), (c) Au nanorod, on a carbon grid, imaged 
with a transmission-electron microscope (TEM). 



30% H 2 2 and 98% H 2 S0 4 , and then coated with an 
MPTMS (3-mercaptopropyltrimethylsilanc) monolayer, 
using a two-step gas-phase silinazation procedure^ The 
functionalized glass substrate is dipped into the Au- 
nanorod solution for 30 minutes. The sample is then 
washed thoroughly in deionized water and dried in air. 
Atomic-force microscopy (AFM) is used to verify that the 
rods are well isolated on the surface, so that individual 
rods can subsequently be probed optically. 

Optical measurements on the single rods are made 
using total-internal-reflection microscopyi2ii22i2& Incident 
light is focussed onto the sample through a glass prism. 
Scattered light is collected using a microscope objective 
and is imaged onto a multimode optical fiber, which se- 
lects a 1.5 /xm spot on the sample for observation. For 
spectral measurements, the light is sent to a spectrome- 
ter equipped with a cooled CCD array detector ( Andor) ; 
for time-resolved measurements, the light is sent to an 
avalanche photodiode (Hamamatsu). 

III. IDENTIFICATION OF SINGLE RODS 

Making single-rod measurements requires carefully en- 
suring that only one rod at a time is being probed. The 
primary method of identifying single rods is to excite 
with incoherent white light and measure the scattering 
spectrum; a typical spectrum is shown in Fig. ^a). The 
narrow resonance, much less broad than the ensemble 
peak, indicates that the scattering comes from a single 
rod. If two or more nanorods are probed, they will have 



different shapes, and thus different plasmon resonance 
frequencies; the measured scattering spectrum will then 
be broader than the single-rod spectrum. 

The scattering spectrum is quantitatively compared to 
a calculation in the quasi-static approximation;2ii22*22i as 
shown in Fig. ^a). We treat the rod as a prolate el- 
lipsoid, and approximate the asymmetric environment of 
the rod as a homogeneous, transparent medium with di- 
electric constant e m = 1.3. For incident light polarized 
parallel to the long axis of the rod, the polarizability of 
the particle is 



e m + L(e - e m ) 

where V is the volume of the rod, e is the dielectric func- 
tion of Au, and L is a geometric factor: 

where e is the eccentricity of the ellipsoid. The scattering 
cross-section is then given by 

k A 2 

Vscat = \a\ , (3) 

where k is the wavenumber of the incident light. 

The imaginary part of the dielectric function of Au is 
taken to bo 27 i 28 . 29 

e 2 (u,T e ) = rf^Ln +4- C (^T e ) , (4) 

w [w$ + i{T e y\ 
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FIG. 2: (Color online) Single-rod scattering signal as a func- 
tion of delay between two incident laser pulses. Left-hand 
side: scattering intensity for overlapping pulses with an en- 
ergy of 47 pJ (light line); envelope of the interference pattern 
(solid red line) ; the same envelope, inverted about the average 
scattering signal at a delay of 75 fs (dashed red line) ; and the 
average of the upper and lower envelopes (heavy line) . Right- 
hand side: scattering intensity for non-overlapping pulses 
with an energy of 94 pJ. Note that both vertical and hori- 
zontal scales are different on the two sides of the graph. 



where u> is the optical frequency and T e is the temper- 
ature of the conduction electrons in the rod. The first 
term is the Drude free-electron contribution; uj p is the 
bulk plasmon frequency, and 7 is the plasmon damping 
rate. The second term is the contribution of transitions 
between the d and the conduction bands. The real part of 
the dielectric function is calculated from this imaginary 
part using the Kramers-Kronig relation. The matrix el- 
ements of the interband transitions and the Drude plas- 
mon frequency are adjusted to reproduce experimental 
dielectric functions^ 

In comparing the calculated and measured scattering 
spectra, the only free parameter is the aspect ratio of the 
rod. For the particular rod in Fig. [Ja), the fitted aspect 
ratio is 5.25, consistent with the rod shapes measured 
by TEM. A different choice of refractive index for the 
surrounding medium changes the fitted nanorod aspect 
ratio, but has no other appreciable effect on calculated 
optical properties. The very good agreement between the 
calculated and measured scattering linewidths is thus a 
clear indication that only a single rod is being probed. 

Identification of single rods is further supported by 
the strong polarization dependence of the scattering, as 
shown in Fig. ^b). The nearly complete modulation is 
consistent with scattering from a single, oriented dipole, 
rather than multiple particles. 



FIG. 3: (Color online) Measured asymmetry in interference 
patterns for three different rods. 



IV. ULTRAFAST NONLINEARITIES 

Nonlinearities of the single nanorods are measured us- 
ing an interferometric scattering technique. The rods 
are excited with 20-fs pulses from a mode-locked, cavity- 
dumped TkSapphire laser^i The pulses are split into 
two equal- intensity parts and the delay of one of 
the pulses is controlled relative to the other by moving 
a retroreflector, using either a calibrated stepper motor 
or a piezoelectric transducer. A single lens focuses the 
two pulses to a common 20-/zm spot on the sample. The 
signal is processed by a lock-in amplifier, which is syn- 
chronized to a chopper that modulates both laser beams. 

Fig. [5] shows an example of the scattering signal from 
a single rod for short delays; the measured interference 
pattern exhibits a pronounced asymmetry in intensity. 
When the laser pulses interfere constructively, the inci- 
dent intensity is doubled, but the amount of scattering 
from the rod increases by less than a factor of two, mean- 
ing that the scattering cross-section of the rod is smaller 
for the higher intensity. The interference patterns do not 
change as the repetition rate of the laser is varied, indi- 
cating that slow, cumulative effects are not important. 
The nonlinearity thus arises within the 20-fs pulse dura- 
tion. 

These measurements make it possible to establish the 
magnitude of the ultrafast nonlinearity. Figure [3] shows 
for different rods the measured asymmetry as a function 
of the laser intensity, /, in each pulse. The nearly linear 
dependence indicates a third-order nonlinearity. That is, 
the scattering cross-section can be written a (I) = + 
la^ , so that the measured asymmetry should be 

^ =4f g (°)-L(3) - ^ 
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FIG. 4: (Color online) (a) Measured single-rod scattering 
spectra for different incident intensities. The rightmost curve 
is the linear spectrum measured using incoherent, broadband 
excitation. The other curves are measured using a single laser 
beam, and correspond to pulse energies of 52, 82, and 192 pj, 
from right to left. The spectra are normalized for ease of 
comparison, (b) Corresponding calculated scattering spectra, 
assuming spectral changes are due to instantaneous heating 
of conduction electrons. 



As shown in Fig. |3 this formula fits our data well. 

The fit to the data gives a normalized nonlinear cross- 
section ct( 3 Vct (0) = (7.5±0.9) x 1CT 11 cm 2 /W. Transient- 
transmission measurements on nanorod ensembles in so- 
lution show a corresponding average nonlinear cross- 
section of approximately 2.5 x 10 -13 cm 2 /W)S the sig- 
nificantly lower value is the result of the large number of 
non-resonant nanorods, whose nonlinearities are smaller 
than and can even oppose those of the resonant rods. 
The single-rod nonlinear cross-section also implies a non- 
linear susceptibility for Au over the laser bandwidth of 
~ 5 x 10~ 18 m 2 /V 2 . As a result of this large non- 
linearity, the change in scattering cross-section can reach 
over 20% at high laser intensities. If the laser power is 
increased further, optical damage occurs, and the scat- 
tering signal gradually and irreversibly decreases. 

Further insight into the measured nonlinearity is ob- 
tained by measuring the dependence of the nanorod scat- 
tering spectrum on incident laser power. Results for a 
particular nanorod are shown in Fig. EJa); an intensity- 
dependent red shift Alo and line broadening A7 are 
clearly seen. Both effects are linear in /; for this rod, 
Aui = 59 meV/nJ, and A7 = 87 meV/nJ. 



V. PICOSECOND NONLINEARITIES 

Having established the magnitude of the nonlinearity, 
we next investigate its time dependence. To do so, we 
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FIG. 5: (Color online) Measured single-rod scattering signal 
as a function of time delay between two incident laser pulses, 
normalized by the measured signal at a delay of 20 ps (points) , 
and calculated change in scattering (lines) . The three curves, 
from top to bottom, correspond to pulse energies of 26, 53, 
and 94 pj. 



perform measurements with longer time delays, so that 
the two laser pulses no longer overlap. Figure [5] shows 
a representative result, and Fig. Ogives similar results 
for different laser powers. The response is characteristic 
of the heating of conduction electrons by the laser pulse, 
followed by their cooling and equilibration with lattice 
phonons Ai2iiLS Increasing the delay up to 150 ps results 
in no detectable change in the scattering signal, indicat- 
ing that effects related to the heating of lattice phonons 
are unimportant on experimental time scales. The data 
in Fig. EJhave been normalized to the signal at long de- 
lays; any ultrafast nonlinearity arising within the laser 
pulse duration will result in a change in this reference 
level, but will not otherwise not affect the time-delay- 
dependent data. 

The picosecond-scale results can be modeled as follows. 
The amount of light transferred from the first laser pulse 
to the conduction electrons is determined by calculating 
the nanorod absorption cross-section; this energy trans- 
fer results in a higher electron temperature, T e . The 
subsequent evolution of T e is calculated by treating the 
conduction electrons and the lattice phonons as two cou- 
pled thermal reservoirs The time evolution of the 
reservoir temperatures is given by 

C e T' e {t) = g[Tt(t)-T e (t)] (6) 
dTl(t) = ff [T e (t)-7Hi)] , (7) 

where C e and Ci are the heat capacities of the electrons 
and the lattice, respectively; TJ is the lattice temperature; 
and g is the electron-phonon coupling coefficient. 
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The effect of elevated T e is to alter the dielectric func- 
tion for Au, as described by Eqn. Q. 27,28 i 29 From the 
dielectric function, the scattering spectrum can be calcu- 
lated according to Eqns. and ©. For higher electron 
temperatures, the plasmon is broadened and red-shifted. 
The modified plasmon response is used to calculate the 
amount of light scattered from the second laser pulse. 

This calculated scattering signal is fit to the data us- 
ing a single free parameter, relating the measured laser 
power to the optical intensity incident on the rod. Fig. 
shows, for a single rod, results for three different laser 
powers. Equally good agreement was obtained for sev- 
eral other laser powers and for other rods, indicating that 
electron heating can account for the measured nonlinear- 
ity on picosecond time scales. 



VI. PLASMON SATURATION 

Unexpectedly, the same thermal model that explains 
the picosecond measurements also quantitatively ex- 
plains the measured nonlinearities on femtosecond time 
scales. More precisely, we can extrapolate the measured 
thermal nonlinearity for a given laser intensity, /, to zero 
time delay; this gives a "dip" D(I) in the normalized 
scattering signal. (See Fig. in this case D ss 9%.) The 
measured values of D(I) can be compared to the asym- 
metries, A(I), of the measured interference patterns. As- 
suming that the only nonlinearity, even for the short- 
est time delays, is due to electron heating, we obtain 
A(I) = 2D (21). (We note that this relation takes into 
account the dependence of the reference level for D(I) on 
the intensity of the two laser pulses.) We observe exactly 
this relation, within our experimental error, meaning that 
we sec no change in the dynamics of the nonlinearity as 
we move from femtosecond to picosecond time scales. 

This observation is consistent with the increase in the 
plasmon linewidth at high pulse energies, as shown in 
Fig. E2a). Fig. Etb) shows theoretical scattering spectra 
for the same pump powers, assuming that the only spec- 
tral changes are due to instantaneous heating of the con- 
duction electrons. Differences between the calculations 
and the measurements, particularly in the line shapes, 
are likely due to femtosecond-scale dynamics that are 
not captured by the extreme assumption of instantaneous 
electron heating. Nonetheless, the good agreement pro- 
vides another strong indication that a nearly thermal dis- 
tribution of electrons is produced in the rod within a time 
short compared to the 20-fs laser pulse duration. 

This is unexpected, since the resonant laser pulses 
should excite plasmons that remain coherent for 15 fs, 
based on the measured linear-scattering linewidth (see 
Fig. Ola)). Any initial nonlinearity would then be due to 
deviation of the plasmon motion from perfect harmonic 
oscillation, caused, for example, by confinement of elec- 
trons by the boundaries of the rod^* The amplitude of 
electron oscillation can be estimated by considering the 



dipole moment D induced by the applied field: 

D = aE = nde , (8) 

where E is the applied field, n is the number of elec- 
trons in the rod, d is their displacement, and e is the 
electronic charge. Using known material parameters for 
Au then gives, for pulse energies of 100 pJ, an electron 
displacement d approximately 8% of the rod length. The 
deviation from simple harmonic oscillation should then 
be considerable, implying a significant coherent nonlin- 
earity. Such a mechanism would also be responsible for 
the generation of third-harmonic radiation by resonantly- 
driven plasmons^*^ 

The absence of any measurable coherent nonlinearity, 
and the immediate emergence of an incoherent thermal 
nonlinearity, indicate that the plasmon cannot be coher- 
ently oscillating over the duration of the laser pulse. In 
other words, the strong, resonant laser excitation must 
be responsible for increasing the plasmon damping rate 
and destroying its coherence. The reduction in the plas- 
mon lifetime then means that it is impossible to resolve 
any coherent nonlinearities with the 20-fs pulses used. 

The increased damping may be due to a greater rate 
of dephasing collisions with the nanorod boundaries or 
to higher-order plasmon-plasmon or plasmon-electron 
interactions. Such interactions are also manifest in the 
significant electron energies that are observed in pho- 
toemission measurements when plasmons are resonantly 
excitedi 36 ! 37 



VII. CONCLUSIONS 

The measurements described in this paper have estab- 
lished for the first time the magnitude of resonant opti- 
cal nonlinearities in single Au nanorods. A benchmark 
value of 20% has been obtained for the nonlinear change 
in the scattering cross-section, using pulse energies that 
induce no optical damage. Surprisingly, resonant exci- 
tation of plasmons results in the same nonlinearity as 
incoherent excitation of conduction electrons. This in- 
dicates that strongly driven plasmons experience a new, 
intensity-dependent damping. There are still several po- 
tential routes towards achieving stronger nonlinearities in 
these systems, such as embedding the nanorods in a po- 
larizablc medium, or assembling them into ordered struc- 
tures. Our current observations thus represent a first step 
towards achieving very large optical nonlinearities on the 
nanometer scale. 
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